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ABSTRACT 


Viscous drag in the transonic regime over an axi-symmetric body with a 
unique aft contour surface is avesipued The forebody is composed of an 
arbitrary ellipsoid. The unique aft contour surface has been obtained by an exact 
solution of the small perturbation transonic equation, using guidelines and tools 
developed at the Naval Postgraduate School. This unique contour allows the delay 
of shock formation in the aft portion, hence delaying the onset of wave drag which 
results in a reduction of the overall transonic pressure drag on the body. The drag 
coefficient thus computed is compared with another axi-symmetric body with the 
same ellipsoid forebody but a simple boat-tailed conical afterbody. Computational 
Fluid Dynamics (CFD) has been used to compute the viscous flow over the two 
bodies using a Navier-Stokes flow-solver. Results obtained confirm the advantage 
of the special shaped afterbody over the conical afterbody by showing the delayed 
formation of shock waves at the aft portion in transonic flow, consequently 
achieving a lower maximum drag coefficient of approximately 5.5%. These results 
can be used in the design of low pressure-drag surfaces for shapes such as 
missiles, projectiles, aircraft external ferry tanks and aircraft engine nacelles for 


improved performance within the transonic flight regime. 
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I. INTRODUCTION 


The study of the aerodynamic characteristics of a 
missile-like object from launch to nominal flight speed will 
typically range from low subsonic to supersonic speeds. 
Entrenched within this range is a flight regime known as the 
transonic flight regime. 

The study of transonic aerodynamics deserves individual 
treatment within the framework of gasdynamics. A typical 
transonic flow with a subsonic free-stream Mach number 
contains a supersonic zone, bounded by the sonic line and a 
shock wave, through which the flow decelerates back to 
subsonic flow. Transonic flows occurs typically within the 
free-stream Mach number range of 0.8 < M, < 1.2, depending on 
the shape of the object [Ref. 1]. 

Transonic aerodynamics is difficult to analyze or predict 
because it must be described by nonlinear equations. The 
transonic flow regime is also complicated by the formation of 
shock waves, the location of which depend on the complex 
interaction of free-stream Mach number and geometry of the 
body. This is in contrast to the purely subsonic or supersonic 
flow regimes where an adequate prediction of aerodynamic 
characteristics can be obtained using linear’ theory. 
Consequently, transonic flows are very sensitive to small 
perturbations in the various flow parameters such as Mach 
number and this makes accurate experimentation difficult [Ref. 
ii 

One of the main aerodynamic forces of interest is drag. 
Objects that achieve flight are designed with aerodynamic 
shapes. The prime objective of such designs is to obtain high 
lift and low drag. A condition associated with transonic flow 
is the significant increase in drag [Ref. 2]. Efforts to 
reduce the drag coefficient in the transonic regime must 
concentrate on reducing the pressure drag (inclusive of wave 








drag) contributions to total drag. During wing design, one 
method to achieve lower drag is to use supercritical airfoil 
sections which are designed to delay the formation of a shock 
wave on the upper surface of the airfoil, thereby enabling a 
higher region of lift generation and a weaker shock to occur 
further downstream of the airfoil upper surface [Ref. 2]. 

A traditional approach to the investigation of 
aerodynamic characteristics is based on wind tunnel test data 
and actual flight test results. Unfortunately, both wind 
tunnel and flight testing are considerably expensive and time 
consuming. Furthermore, wind tunnel testing in the transonic 
regime requires special wind tunnels with adaptive walls in 
order to nullify possible reflected shock-interaction between 
the tunnel walls and the test object within the test section 
[Ref. 1]. As mentioned in Reference 1, it is also difficult to 
obtain highly reliable data from wind tunnel experimentations 
in the transonic flow regime. 

In contrast, there has been much progress recently in the 
field of Computational Fluid Dynamics (CFD). The judicious use 
of CFD, in combination with the appropriate modelling, has 
- enabled aerodynamicists to predict aerodynamic characteristics 
within a faster time frame and at relatively lower costs. 
Numerous codes are now available that run as both the Euler 
and Navier-Stokes flow-solvers [Ref. 3]. Many institutions are 
now using CFD as a "first-cut" predictor for improvements in 
aerodynamic designs Sy a6: “tol to extrapolate aerodynamic 
characteristics from known design performances. With the 
increased experience in CFD, such as better turbulence 
modelling and more efficient and robust codes, and with modern 
developments in high speed computing, one can expect to see an 
even larger proportion of aerodynamic design being 


accomplished using CFD. 





II. TRANSONIC FLOW 


Transonic flows are characterized by the simultaneous 
presence within the flow field of both subsonic and supersonic 
regions. The properties of transonic flows can be described 
using the equations of gas dynamics, namely, the equations of 
state, continuity, momentum and energy. If the object of 
interest is a slender body, we can use small perturbation 
theory to simplify the flow over the body [Ref. 2]. In fact, 
small perturbation theory can be applied throughout the entire 
Mach number range from low subsonic to supersonic speeds. 
Typically, the derivation is based on assumptions that the 
flow is steady, irrotational and isentropic with no energy 
transfer, no body forces and no shear stresses (1.e., inviscid 
flow). 


A. SMALL PERTURBATION THEORY 


For steady, isentropic and inviscid flow, the governing 


equations are 


i) Continuity Equation 
O(pu) , O(pv) , d(pw) _, (1) 
Ox dy OZ 


Expanding and re-arranging equation (1) 


Ou, Ov Ow __ op Op, op 
Pa: Pap ce tox on en. (2) 


ii) The Euler momentum equations are given by 


av, ,av,,,8v__ 1 8p 
Uax dy oz p oy 
aw, dw, __ 2 ap 
Ox “Sy "Oz 9 OZ 





(3) 


(4) 


(5) 


Since the flow is isentropic, we can express the speed of 


sound as 


(6) 


Multiplying the x-momentum equation (Eqn 3) by u yields 


u? $2 suv suwle .._Uu Op Op 


oy Oz p Ox dp 


Substituting Egn (6) into Eqn (7) yields 


u2 Se suv ae Spy len IEP 2 


Oy OZ p Ox 


Similarly, by multiplying the y-momentum and 


equations by v and w respectively yields 


OV eat OV gy Vig OR a2 


Ox Oy Oz p oy 


uWSy Oy « OZ~COp Oz 


(7) 


(8) 


z-momentum 


(9) 


(10) 





re 


Summing up Eqns (8), (9) and (10) gives 


2 BU gy BU arr DY ayy OY ay? Ya yy DY ow DY + vw OY +2 OH = 
‘ 5x OV Gy ae OV ax a Oe on Oz 
aja® (uP. +y Ph +w2l) (11) 
a) ox Oy az 
Substituting Eqn (2) into Eqn (11) yields 
2 Ou Ou Qu OV 42 OV 4 OV 4 ry OW + vw OW + > OW _ 


Wo ay Ma Max” Oy Oz.)Ox)SCO yO 


Ou, Ov, Ow (12) 


By rearranging Eqn (12), we obtain 


du_u? du), dv_v? dv), ( aw_w? dw) 
Ox a2 Ox Oy a2 oy OZ g2 02 


_ (uv du, uv Ov) _ (uw du, uw dw 
a? Ody a? Ox a2? OZ a2 Ox 
(YW OV, Vw OW). 
ede ae ay" (13) 
Since the flow is irrotational, we have 
Ou_ ov 
frdiceay odd 4 
ay Ox (14) 
du _ Ow 
ao ox (15) 
Ow _ OV 
By Oa (16) 








Substituting Eqns (14), (15) and (16) into Eqn (13) 


rearranging yields 


Ou ,,_u*,,0V/,_v*,, Ow. _ Ww? 
By it a2 Oy toe a2? 


a? Oy a2 Ox = a2 Oz 


For a slender body moving in the x-direction, 
perturbation theory gives the velocity field as 


u=U.tu’ ; v=v’ ; w=w’ 


Since 


u>u’ ; urv’; usw! 


We can simplify Eqn (17) to 


_u*, du, dv, dw_ 


Se eee — 


a’ Ox Oy OZ 





en 


and 


(17) 


small 


(18) 


(19) 


(20) 


iii) The flow is isentropic and the energy equation can be 


expressed as 


2 
H,=h+ — =Constant 


We have 


VT V 
6 es ee eee 
~ 2 2 


We can use the isentropic relationships 


=—YR_ ; a?=yrt; m=+ 
(Y>L) a 


(21) 


(22) 


(23) 


Substituting Eqns (23) into Eqn (22) and rearranging yields 


a?=a2-1=* (V2-vVé) 





2_tr2 
aii) ( “a Vo ) m2 
a2 2 Ve 


Since the flow is 3 dimensional, we have 
V2=u2t+vet+w? ; Vi=U. 

Substituting in Eqns (18) and rearranging yields 
V2-U2=2U,u!tu*+v"?+w"? 

Substituting Eqn (27) into Eqn (25) yields 


1. 2U.u'tu@+v@+w” 
fies per Ok ee a 


: é 
2 2 2 ) Ma 
Ae Us 
Taking the reciprocal of Eqn (28) 
ae al 
/ f2 f2 f2 
Ae fis gece ag © Sila oh di ) me 
2 8 


Using the Binomial Expansion on Eqn (29) 


2 = fio Pi he 
Ge 2 Mol yy? | 2Ueu tu tw) 4H. O.T. 
a? 2 U2 


From Eqn (20), the first term can be expanded as 





(1- u- ae U2+2U,tu” Uza? =1-( U2+2U.+u? , az 
a2 a? Vaz ae Mee ee 


(24) 


(25) 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 








Substituting Eqn (30) into Eqn (31) yields 








20 ee Vol yey cu Uv 
—)=1-(1+ Pees MO tL + + 
ee -) ( ia <i [ 5 Mz ( a ae oe 


f2 f2 f2 
,W 


8 tw 


u 


Eqn (32) can be simplified by neglecting the higher order 


terms (H.O.T.) to give 


(1- a =1-M2 [1+ 24 u(i+Tm)) 
a? UU, 


Substituting Eqn (33) into Eqn (20) yields 


14-22 y 1 du, Ov, dw 
[1 -Me- SMe (14-5 Me)] at ay az ° 


Rearranging Eqn (34) yields 


du. Ov, dw_ 2u’ Pie a Cu 
(1 “Ma) + By Oz. Spe (1+ 5 Me) a 


(33) 


(34) 


(35) 


Rewriting Eqn (35) in terms of the perturbation velocities 


gives 
O{U,+u") | av’, dw’_ 2u’ -1 d(U,+u") 
SNS Vo a eh ins © oy PRON AEN 


We assume that 


(37) 


Expanding Eqn (36), substituting into Eqn (37) and neglecting 


2nd order perturbation terms yields 


_ du’, Ov’, Ow’ _ M. Y=2 y2) 2y/9U! 
(a8) Bs Be AE (20 seal le 3 


(38) 





Since the flow is irrotational, we have 


1 OD . 1. DD . y/- D 
U ox) © oy’ ” Oz 32 


Substituting Eqn (39) into Eqn (38) gives’ the small 
perturbation potential equation for transonic flows 


2 
(1-M2) bye yy Dag" (14 GME) 26,9 (40) 
U2 2 


B. AN EXACT SOLUTION FOR AXI-SYMMETRIC BODIES IN TRANSONIC 
FLOW 

Solutions to the two-dimensional, small perturbation 
potential equation for transonic flow have been accomplished 
by Biblarz [Ref. 4] by using the separation of variables 
approach. The solution of the small perturbation transonic 
flow equation has enabled the generation of shockless boundary 
surfaces. Al-Hashel [Ref. 5] reports on axi-symmetric boundary 
surfaces for M=1.0, 1.10 and 1.20. 


LO 








III. DRAG 


Drag can be defined as the fluid dynamic force acting 
against an object placed in the stream of airflow. Basically, 
the drag force can be divided into two broad categories. The 
categories are viscous or skin friction drag and pressure 
drag. 

Viscosity is the molecular "“resistance" which fluid 
particles exhibit against displacement in relation to each 
other and with respect to the surface of solid objects. Such 
resistance presents itself in the form of frictional force, 
which is comparable to that of solid surfaces sliding along 
each other. This frictional force is a tangential force which 
occurs when a fluid (such as air) flows past a body. This 
force is the skin friction drag. At very low speeds or in 
flows associated with fluids of high viscosity and bodies of 
very small dimensions (1.e., at very low Reynolds number 
flow), viscosity is the predominant property determining the 
drag of a body [Ref. 6]. At higher Reynolds number flows, 
Other forms of drag; such as turbulent skin friction drag, 
begins to have significant contribution. 

In contrast to the skin friction drag (which is a 
tangential force), pressure drag results from the distribution 
of pressure forces normal to the body surface. Within the 
subsonic range of speeds, inviscid fluid dynamic theory 
predicts that the flow will close in behind an object, without 
any losses. Positive pressures at the front of the object are 
counterbalanced by the negative pressure of equal magnitude on 
the rear, resulting in the pressure drag being zero. This 
phenomenon 1S known as D’Alembert’s Paradox. In real flows, 
however, such a phenomenon does not exist. The pressure at the 
rear side of an object can be very much lower than predicted 
from non-viscous flow conditions. Viscosity causes the 


formation of a boundary layer which may cause flow separation 


-- 
-- 








on the surface of an object. The resultant pressure 
differential between the pressures on the forward side and the 
lower pressures on the rear side represents the pressure drag. 
In blune or bluff bodies, this type of drag 1s usually much 
higher: than the skin friction drag [Ref. 6]. For streamlined 
or aerodynamically shaped bodies, the reverse is true. For 
flow in the transonic regime, the formation of shocks 
contributes to the pressure drag due to the pressure 
differential between the forward portion of an object and the 
rear portion of the object. Figure (3.1) [{Ref. 7] shows the 
typical drag composition of an axi-symmetric projectile shaped 
body. Figure (3.2) [Ref. 8] shows the variation in skin 
friction drag and wave drag as a function of free-stream Mach 


number for a body of revolution at zero incidence. 


Wave Drag 


CDo 


Base Drag 


Subsonic Pressure Drag 


Skin Friction Drag 





06 0.7 0.8 09 1 4.4 1.2 1.3 1.4 4.5 1.8 
Mach Number 


Figure 3.1. Components of Drag [Ref. 7] 
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Friction Drag 


Wave Drag 


Coefficient of Drag 





Figure 3.2. Friction and Wave Drag Variation [Ref. 8] 


One well-known method of reducing the drag on an axi- 
symmetric body is by streamlining the afterbody, i.e., 
reducing its base diameter gradually. Figure (3.3) [Ref. 6] 
shows the reduction in zero incidence drag coefficient for 
boat-tailed afterbodies. Other approaches of reducing drag on 
a streamlined body includes mechanical methods such as vortex 
generators and methods for controlling the boundary layer such 


as the boundary layer suction method. 
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Me= 0.8 ; Re = 6 x 10° 
Direction of Flow ----> 


—}. oo etree fh Coo = 0.17 
—- S- Co, = 0.08 


Figure 3.3. Drag Coefficient of Boat-Tailed Afterbodies [Ref. 
6] 


0.06 


Numerous methods are available to calculate the 
theoretical values of drag on an object. However, there are no 
totally satisfactory method of accurately calculating viscous 
drag, particularly transonic drag. This is due in large to the 
necessary assumptions and simplifications to the full Navier- 
Stokes equations to facilitate mathematical solutions and the 
need to use turbulence models for the boundary layer effects 
in high Reynolds Number flows. In this work, both the skin 
friction drag and the pressure drag can be obtained directly 
from the results of CFD computations [Ref. 5 9 ee 
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IV. COMPUTATIONAL FLUID DYNAMICS (CFD) 


Rapid advancement in the speed and efficiency of 
computing and better modellings has led to the emergence of 
computational fluid dynamics as an important branch of 
aerodynamics. CFD complements the experimental and theoretical 
branches of fluid dynamics by providing a relatively accurate 
and cost-effective means of simulating real flows. Another 
significant advantage of CFD is the convenience of being able 
to vary any relevant parameter of the flow and also the 
ability to "switch off" specific terms in the governing 
equations so as to assist the researcher in understanding 
their contributions to the resultant flow. 

In this work, CFD was used to compute the axi-symmetric 
flow over complete body models which are composed of a 
forebody (ellipsoid) and an afterbody (conical or special 
shape) using a Navier-Stokes flow-solver. The surface profile 
of the special shape afterbody was reported on by Al-Hashel 
[Ref. > using an exact solution method for small 
perturbation. The conical afterbody has a base diameter ratio 
fa. dG.) Of 0275, and-@ conical. turning angle (fp) of 14 degrees: 
A generic forebody in the shape of an ellipsoid was patched 
onto the afterbodies to establish the CFD working models. 


A. GRID GENERATION 


A computational grid system is a necessary part of any 
numerical solution based on a finite-difference or finite- 
element method. The selection of a grid system is based 
primarily on the requirement for accuracy in the final 
solution. Secondary considerations are the effects on 
computational efficiency of the solution algorithm using 
available computer architecture and, finally, the ease of grid 


generation. 


io 





The grids for the two models used were generated by 
Priyono [Ref. 9] using the computer programs GRAPE [Ref. 10] 
and GRIDGEN2D [Ref. 11]. These programs were used to generate 
two-dimensional grids which were subsequently converted into 
3-D grids utilizing the FORTRAN codes d2d3.f and rotategr.f 
[Appendix B]. GRAPE was used to generate the grid for the 
model with the special shape afterbody (model CNTR11). 
GRIDGEN2D was used to generate the grid for the model with the 
conical afterbody (model CBXCN2), Since GRAPE cannot 
accomodate the non-smooth points associated with the conical 
afterbody. The grids for the two models are shown at Figures 
(1) and (2) in Appendix A. In this work, the correspondence 
between the two body profiles was verified priror to turener 
analysis. This correspondence is shown in Figure (3) in 
Appendix A. 

Since the analysis was for the complete body, an O-type 
grid was used. Previous work by Priyono demonstrated the need 
to use a fine grid due to the significant propagation of the 
bow shock above and below the body in transonic flow, 
especially in the vicinity of freestream Mach number of unity. 
It was necessary for the bow shock to diminish before it 
reached the grid boundary, thereby ensuring that the shock 1s 
not reflected by the grid boundary and contaminated the 
solution. The outer boundary of the fine grid was 5 times the 
body's length. The grids were generated only for half of the 
complete body taking advantage of symmetry. The grid size used 
ic 7 7RLSR LO. 


B. FLOW-SOLVER 


The OVERFLOW program was developed by NASA Ames Research 
Center. The version used is Version 1.6ag released on 30 April 
1993 [Ref. 3]. The program uses either the 3-D Euler or 
Navier-Stokes flow-solvers for inviscid or viscous flow 


computations respectively. The selection of Euler or Navier- 
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Stokes flow-solver is accomplished by setting the parameter 
VISINP (viscosity input) to True or False in the input file. 

Before applying the flow-solver code (OVERFLOW), the 
formatted 3-D grid file (grid.for) must be converted into an 
unformatted input grid file (grid.in). This is accomplished 
using the FORTRAN code readx.f [Appendix B]. An input file 
(overflow. in) containing all the input parameters must then be 
written for running OVERFLOW. The input parameters consist of 
the number of iterations, timesteps, calculation methods, 
smoothing method, types of in-flow and out-flow, boundary 
conditions and types of turbulence models. The value of angle 
of attack (ALPHA) depends on the orientation of the grid in 
the coordinate system. In our case, ALPHA was set at 180° 
(i1.e., flow comes from the x-positive to the x-negative 
direction). A typical input file (overflow.in) is provided in 
Appendix B. 

In his work, Priyono investigated the flow past 
afterbodies only using both the Fuler and Navier-Stokes flow- 
solvers. Priyono also investigated the flow past complete 
bodies but only with the Euler flow-solver. For the present 
work which investigates the flow past complete bodies using 
the Navier-Stokes flow-solver, the effects of viscosity have 
to be “switched on". This 1s accomplished by setting VISCJ, 
VISCK and VISCL parameters to TRUE. This in effect includes 
viscous terms into the J (along the body in the x-direction), 
K (around the body or x-axis) and L (outwards from the body 
towards the outer grid boundary) computational directions. The 
number of turbulent wall regions NTURB is set to 2. This 
corresponds to the turbulence models ITTYP of 1 and. 11, 22664, 
the Baldwin-Lomax boundary layer model (with variable Degani- 
Schiff cutoff) for the body in general and the Baldwin-Lomax 
shear layer model (with variable C,,,) to model the wake behind 
the body. The basis of selecting the turbulence models comes 


from the OVERFLOW manual which recommends the use of the 
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Baldwin-Lomax model for flows with mild separation. 

The OVERFLOW program generates output files such as 
overflow.out, fomo.out, rpmin.out, resid.out and q.save. The 
file overflow.out provides information on the input parameters 
which are used in the computational run. The file fomo.out 
provides the force and moment history, which in this case, 
includes coefficients of drag due to pressure (which was also 
computed when using the Euler flow solver) and viscosity 
(i.e., skin friction drag). The file rpmin.out provides the 
minimum density/pressure/ratio of specific heats history. The 
file resid.out provides flow solver residual history for the 
five primary variables (i.e., p, pu, pv, pw and e, where e 1s 
the total energy). The variable that is used to_ show 
convergence of the solution is the density residual. In 
general, the convergence criterion can be defined as the 
reduction in the density residuals by two orders of magnitude 
and the trend towards further reduction. The file q.save is 
the solution file in PLOT3D format [Ref. 12] when OVERFLOW has 
run satisfactorily. If the flow-solver detects a case of 
negative density or pressure during the computation process, 
an output file called q.bomb will be generated instead of the 
solution output file q.save. The file q.bomb thus represents 
the solution prior to the error being detected. 

The output from resid.out is a text file and the result 
can be plotted using a plotting software such as GNUPLOT 
(Ref. 13]. The output from fomo.out is also a text file and 
the values of coefficient of drag due to pressure and 
viscosity can be read off directly at the end of the file. The 
coefficient of drag can also be plotted using GNUPLOT to show 
convergence to a particular value. The output from q.save 1s 
an unformated file which needs to be converted into a 
formatted file (q.form) using the FORTRAN code readq.f 
[Appendix B]. The formatted solution file q.form can be 
plotted using either PLOT3D or FAST [Ref. 14]. Both software 
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packages are capable of generating plots showing the Mach 
number distribution, shock locations due to pressure 
difference, velocity vector distribution and many other 
options of flow visualization over the bodies. In this work, 
PLOT3D was used for all the flow visualization plots. 


C. CFD SOLUTIONS 


Initial turbulent flow computations were accomplished 
with the setting of ITDIR of 3 andl, corresponding to the L 
and J coordinate directions away from the wall or shear layer 
for the Baldwin-Lomax boundary layer model and the Baldwin- 
Lomax shear layer model respectively. Solutions obtained from 
the computations displayed an asymmetry at the wake of the 
bodies. Logically, there should be no asymmetry Since the 
computations were for zero incidence flow. 

Re-examination of the input for ITDIR strongly indicated 
that the parameter should be set to 3 and 3, corresponding to 
the coordinate direction of L for both the turbulence models. 
Subsequent computations with this setting produced symmetry at 
the wake. The output also displayed a clear momentum deficit 
within the wake region aft of the bodies. 

Having resolved the parameters for the turbulence 
modelling, the next step was to ensure effective convergence 
of the solutions. A slow start was used with DT = 0.01 and 
CFLMIN = 3.0 for all the runs. After 2000 iterations, it was 
observed that the density residuals were only converging to 
approximately 1.5 orders of magnitude. Since the viscous 
solution includes skin friction drag, the residual would have 
to show at least 2 orders of magnitude of convergence. In 
order to increase the convergence rate, the runs were 
initiated for either 500 or 1000 iterations with a CFLMIN = 3, 
depending on the free-stream Mach number. Subsequently, 
restarts were used and the values of CFLMIN or DT were varied. 
Typical variations of CFLMIN and DT are tabulated at Table (1) 
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in Appendix C. Figures (4) and (5) in Appendix A show typical 
plots for density residuals. The fluctuations in the residuals 
are caused by the effects of turbulence modelling. 

The critical free-stream Mach number for both bodies was 
computed to be approximately past 0.8. This is in contrast to 
the Euler results reported by Priyono which show the inviscid 
critical Mach number as approximately 0.7 and that sonic 
conditions were reached in the regions where the forebody and 
afterbodies mate with the mid-sections. Figures (6) and (7) in 
Appendix A shows the Mach contours for CNTR11 and CBXCN2 
respectively at M=0.8 and the highest local Mach number 
displayed is 0.9, indicating that the critical Mach number has 
not been reached. Figures (8) and (9) in Appendix A shows the 
Mach contours and shock locations for CNTR11 at M,=0.95 
respectively. Comparing Figure (9) with Figure (70). an 
Appendix A, which shows the shock locations for CBXCN2 at 
M.=0.95, note that the rear shock is further back at the aft 
portion of CNTR11. The viscous critical Mach number was higher 
due to the presence of the boundary layer which lessened the 
acceleration of the airflow over the forebody and also had a 
"rounding" effect on the afterbodies which translates into a 
slower acceleration around the aft regions. 

Since this work uses the Navier-Stokes flow-solver, it is 
important to verify that the resultant flow is indeed a 
viscous solution. The plot at Figure (11) in Appendix A shows 
the velocity vector distribution over CNTR1i1 at a free-stream 
Mach number of 1.2. It is important to note the distance 
boundary layer in the plot. It was also necessary to validate 
the no-slip condition on the surface of the bodies. The 
appropriate input parameters were used in PLOT3D to show only 
the Mach contours on the surface of the bodies. The results 
show a Mach contour of zero on the surface of the bodies, 


thereby confirming the no-slip COnadl ELON: 
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The plots showing the shock locations based on pressure 
gradient complement the plots of Mach number distribution. The 
plots at Figures (12) and (13) in Appendix A for CNTRI1 and 
CBXCN2 respectively at M.=1.05 show that a bow shock has 
formed once the free-stream Mach number goes supersonic. It is 
interesting to note that the plot at Figure (10) in Appendix 
A for free-stream Mach number of 0.95 corresponds to a 
schlieren picture at Figure (14) in Appendix A of a similar 
flow over an axi-symmetric body at free-stream Mach number of 
0.95 at zero incidence [Ref. 15]. 
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V. RESULTS, CONCLUSIONS AND RECOMMENDATIONS 


A. RESULTS 


The values of drag coefficient are obtained directly from 
the OVERFLOW output files fomo.out. The output from these 
files 1s plotted to show that the values of pressure drag 
coefficient and skin-friction drag coefficient are converged. 
A typical plot is shown at Figure (5.1). The viscous drag 
coefficients of the two bodies at the various Mach numbers are 
compared in Figure (5.2). The corresponding inviscid drag 
coefficients as computed by Priyono are indicated in Figure 
(5.63) 
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Figure 5.1. Convergence of Drag Coefficient 
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As expected, the viscous drag coefficients are higher 
than the inviscid drag coefficients. These results show that 
skin friction drag coefficients are significant for the 
viscous solutions. The viscous solutions show a maximum drag 
coefficient difference of approximately 5.5% between the two 
bodies. In comparison, the Priyono [Ref. 9] reported a maximum 
drag coefficient difference of approximately 15% between the 
two bodies for the Euler solutions. The likely difference 
between these is the smoothing of the afterbody by the 
boundary layers. 

The pressure drag coefficient remains relatively constant 
from the sub-sonic free-stream Mach number regime up until the 
Critical Mach number. Subsequently, the pressure drag 
increases due to the addition of wave drag resulting from the 
formation of shock waves. The skin friction drag coefficient, 
on the other hand, decreases with increasing free-stream Mach 
number. The drag coefficients as shown in Figure (5.2) is in 
good agreement with typical values as shown in Figures (3.1) 
and (3.2). Preliminary results for free-stream Mach numbers of 
1.3 and 1.5 show that the drag coefficient decreases from the 
peak value shown at Figure (5.2). 

The results show that the special shape afterbody did 
indeed provide a lower drag coefficient as compared to the 
conical afterbody. The plots from shock location for free- 
stream Mach number of 0.95 clearly shows the position of the 
rear shock wave as being further aft of the special shape 
afterbody when compared with the rear shock wave location for 
CBXCN2. 


B. CONCLUSIONS AND RECOMMENDATIONS 


The OVERFLOW code has proven to be a useful CFD tool for 
this work. In hindsight, it would have been preferable to 
generate the grids for both the bodies using one grid 
generation method instead of using GRAPE for CNTR11 and 
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GRIDGEN2D for CBXCN2. The advantage would have been complete 
commonality between the grids for the two bodies with the only 
variations being at the aft portions. As mentioned earlier, 
although the profile of the two bodies were verified, the 
positioning of the grid points on the two bodies (especially 
on the forebody) were different. This difference may affect 
the accuracy of the solutions. 

The investigation of the two axi-symmetric bodies using 
CFD has shown an advantage of the special shaped afterbody a 9 
the form of lower transonic viscous drag coefficient over the 
conical shaped afterbody by a maximum difference OF 
approximately 5.5%. This advantage is over and above the 
benefits that one can achieve by boat-tailing the afterbody as 
shown in Figure (3.3). These results can be used to design low 
pressure drag surfaces for shapes such as missiles, 
projectiles, aircraft external ferry tanks and aircraft engine 
nacelles for improved performance within the transonic flight 
regime. 

Finally, this work may be continued with experimental 
validation of the advantage of the special shaped afterbody 
over the simple conical afterbody using an appropriate 
transonic wind-tunnel. Further CFD work could also be carried 
out to investigate the effects of angle-of-attack on the two 
bodies. Such efforts would contribute useful experimental 
data on the transonic flow regime, the likes of which is 


certainly lacking currently. 
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APPENDIX B. FORTRAN PROGRAMS AND OVERFLOW INPUT FILE 


4> 


Pare aaa a ara are war wae ear ean eee eee ee ee ne ee ee ee eee eee ee 
* 


* 


* 





D2D3 .F 


* This FORTRAN program will prepare a formatted 2-D grid 
* data file for processing by ROTATEGR.F 


Parra ar eran area arena ara arr rar arararar ara aera ee See Gen ee ee ee ee ee ee ee ee 


+ 


+ 
+ 


* 
* 
* 
* 
* 
a 


WRITE (14, *) 


real X(77,120) ,Y(77,120) ,2(77,120) 


open (unit=12, file=’grid2d.for'’,status=’unknown’ ) 
open (unit=14, file=’cone3d. for’, status=‘new’ ) 


READ(12,*) IDIM,JDIM 


READ(12,*) ((X(I,J),I=1,IDIM),J=1,JDIM) , 
((Y(I,J),I=1,IDIM) , J=1,JDIM) 


Z2(I,Jd)=0.0 
WRITE(14,*) IDIM,JDIM,1 


STOP 
END 


ROTATEGR.F 


This FORTRAN program takes a formatted 2-D grid and 
generates an Axi-Symmetric an unformatted 3-D grid 


Converts grid2d.for to cone.for 


dimension x(77,120), VI 204, 


dimension xx(77,15,120), yy(77,15,120), 


character*30 f£ni 

character*30 fno 

print *, ‘Input 2-D grid filename’ 
read (*,21) fni 

rewind 2 

open (2, file=fni, form=' formatted’ ) 
read (2,*) 11,31,k1 


read (2,*) (( x(1,3),1=1,11),3=1,351), 


CC VO oe bey ye os 
( ( 23). PEN ad) eda I) 
close(2) 
pi = 4.*atan(1.) 
print *, ‘No of Planes in j dir ? ’ 


read (5,*) ym 
dth = (180. / (jm-1)) * (pi/180.) 
do 11 1 2+41,11 


do 11 j 1,31 

ko9 

sexi 2.k) = xia) 
yylig2;k) = 04 

Zook) = ela) 
continue 46 


((X(I,J),I=1,IDIM) ,J=1,JDIM), 
((0.0,I=1,IDIM) ,J=1,JDIM), 
((Y¥(I,J),I=1, IDIM) ,J=1,JDIM) 


zz(77,15,120) 


eR Rt RR KR KR RK RRR KKH KKK KKKAHKEHAKAKAAKKKaKaAKaEKKRA KEK 


* 
* 
* 
* 


rr ae ae ee ee ee ae ae ee a ee ee ee 


do 20 j=3,jm+tl1 

do 20 i=1,im 

do 20 k=1,km 
xx(i,j,k) = x(i,k) 
th = (3-2) *dth 


yy(i,j,k) = sin(th)*z(i,k) 
zz(i,j,k) = cos(th)*z(1i,k) 
continue 

do 30 i=1,1im 

do 30 k=1,km 

j=1 

xx(i,j,k) = xx(1i,3+2,k) 
yy(i,j,k) =-yy(1,3+2,k) 
22(3,3,k) = 22(1,3+4+2,k) 
jJ=jym+2 

xx(1,3,k) = xx(1,j-2,k) 
yy(i,j,k) = -yy(1i,j3-2,k) 
zz(i,j,k) = zz(1,j3-2,k) 
continue 

jm = ym+2 


print *, ‘Output 3-D grid filename =' 

read (5,21) fno 

rewind 3 

open ( 3, file=fno, form=’unformatted’ ) 

write (3) im, jm,km 

write (3) ((( xx(i,j,k), i=1l,im), j=1,3jm),k=1,km), 
+ ((( yy(i,j,k), 221,1im), j=1,3m),k=1,km), 
+ ((( zz(i,j,k), i2=1,im), j=1,3m),k=1, km) 
close (3) 

stop 

end 
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Se ete HH TKK IK III IK IH IIIA IIIT KIKI HIRI KIKI KIRKE KEK EEK EKEEKREK KE A 
* 


* READX .F 

7 This FORTRAN program will read formatted data files and * 
* convert them to unformatted binary data files for input 
* into OVERFLOW - Converts grid.for to grid.in * 


eH KKK TR HMI IK III IIIA IIA IA IARI IIIA HAITI KKK KEK EAH ERE EKEEKEKERE 


real x(77,15,120),y(77,15,120) ,2(77,15,120) 

open(unit=12, file='grid.for’,status=’unknown’ } 

open (unit=14, file='grid.in’,status='new’, form=‘unformatted’ ) 
read(12,*) IDIM, JDIM, KDIM 

read(12,*) (((X(I,J,K),I=1,IDIM) ,J=1,JDIM) ,K=1,KDIM), 


+ (((Y¥(I,0,K),I=1,IDIM) ,J=1,JDIM) ,K=1,KDIM), 
+ (((Z(1,3,K),I=1,IDIM) ,J=1,JDIM) , K=1,KDIM) 
write(14) IDIM, JDIM, KDIM 

write(14) (((X(1,3,K),I=1,IDIM) ,J=1,IDIM) , K=1,KDIM), 
+ (((Y(I,d,K),I=1,IDIM) ,J=1,JDIM) ,K=1,KDIM), 
+ (((Z(I,3,K),I=1,IDIM) ,J=1,JDIM) ,K=1,KDIM) 
stop 

end 


eK KK KKK KKK KKK KKK KK KKK IK KK TK TK KKK IK IK KTR KKK KKK KK KKK KK KS 
ia READO.F * 
* This FORTRAN program will read unformatted binary data . 
* files and convert them to formatted data files for input 


* into PLOT3D - Converts q.save to q.form 
Bn KK KK IK KKK KK KKK KKK TK KK KK KKK TRIKE KK KKK KKK KKK KKK RK KK KKH 


* 
* 


dimension q(77,15,120,5) 

open (unit=1, file='q.save’,status=’unknown’, form=‘unformatted’ ) 
open(unit=20,file=’q.form’, status='new’, form=’ formatted’ ) 
read(1' ni, nj, nk 

read(i)fsmach,alpha,re,time 
read(1)((((q(i,j,k,nx),i=1,ni),j=1,nj),k=1,nk),nx=1,5) 
write(20,*) ni, nj, nk 

write(20,*)fsmach, alpha, re, time 
write(20,*)((((q(i,j,k,nx),i=1,ni),j=1,nj),k=1,nk),nx=1,5) 
print*, ni, nj, nk 

stop 

end 
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INPUT FILE FOR OVERFLOW (OVERFLOW. IN) 
KH KKK KKK KAKEKKEKEKKEKAKKKKKKKKKKKS 
SGLOBAL 
CHIMRA 
NOT 
SEND 
SFLOINP 
ALPHA 
SEND 
SVARGAM 
_ SEND 
SGRDNAM 
NAME = ‘Half complete body, CNTR11 , 77x15x120 Grid’, 
SEND 
SNITERS 


a NSTEPS = 1000, RESTRT ee NSAVE = 100, 


wou 
rr] 


0, 


6.Q00E6, TINF 520.0, 


180.0, FSMACH = 0.95, REY 


IRHS = Q, ILHS = 2, IDISS = 2, 


= 0.01, ITIME = 1 TFOSO = 1.00, CFLMIN = 3.0, 


lv, 
a 


O02; 
0.02, 
O02. 


.00, DIS4J 
DIS4K 
.00, DIS4L 


WoW ft 


H 
NY 
Ae | 
ty 
O 
iy 
Hou Ww Wa 
OrRNN ND 
Oo 
a 
ep) 
tO 
x 
Wott ot 
tO b&b 
>) 
© 


SEND 
SVISINP 
VISCJI 
NTURB 
LTTYS 
TrDIR 
JiisS 
JTLE 
KTILS 
KTLE 
LTLS 
LTLE 
TLPARI 
SEND 
SBCINP 
NBC 
TBI? 
IBDIR 
JBCS 
JBCE 
KBCS 
KBCE 
LBCS 
LBCE 
SEND 
SSCEINP 
SEND 49 


me VIiSCK @ 2225  ViSChoeS <P, 


Yow Ww RM 
| 


. be 
WRPrRPUPRPERPHPWHDND.- 


~ ~ ~ ™~ ™ ~~ ~~ ~ ™~ ~ ~ 
= 


foi i ol 
{ 


© 


ey 


HUE RRR UO 
~J 
~J 


pat 
yn 
— 
a" 
uy 
ay 


4 
~] 
~J 
~J 
~J 
~] 
~J 
~] 


a 
a 
bell 
Me UTR JR WU 


ay 
NO 
© 


120, 120, 120, 120, 
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APPENDIX C. TABLE OF OVERFLOW INPUTS 


oul 











ep) 
mm 
O 
rah 
C4 
ee 
[x] 
HH 
ti 


Ud 


0 


CBXCN2 





OL 


0 


7.0 500 


1000 


<U 


3 


Ol 


0 


we 


1 


CBXCN2 


200 


300 


500 


1000 


Ol 


CNTR 


O01 


QO. 


Pa @ 0 1000 


so 


0 


CNTPRID 


800 


me 


5 


<2 © 
© a 
N GN 
© © 
N ~ 


co a 
© © 


0.005 


1000 
500 


0 


O01 3 


CEs, 


al 


CNTR11 


TABLE 1 
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